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We study the algebra of constraints of quantum gravity in the loop representation based on 
■ Ashtekar's new variables. We show by direct computation that the quantum commutator algebra 

reproduces the classical Poisson bracket one, in the limit in which regulators are removed. The 
calculation illustrates the use of several computational techniques for the loop representation. 
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An important element in the classical formulation of a canonical theory which has constraints is their algebra. If 
taking Poisson brackets of the constraints of the theory leads to quantities that can be expressed as combinations 
of the constraints, these are said to be first class. If the constraints are first class, quantization can proceed by 
O" 1 requiring that the quantum operator version of the constraints annihilate the wavefunctions. Upon quantization, the 
^ | Poisson bracket algebra of constraints should translate into an algebra of commutators where the commutator of any 
constraints should again be expressible as a combination of the constraints. The combination may involve coefficients 
| that are functions of the canonical variables (which become operators at a quantum level). It is therefore crucial that 
such operators appear to the left. If this were not the case, imposing that the wavefunctions be annihilated by the 
constraints could imply — via commutation relations — extra conditions on the wavefunctions. The quantum theory 
found would therefore not necessarily correspond to the classical theory one started from. 

From a physical point of view, constraints generate via canonical transformations the symmetries of the theory. 
Therefore respecting the constraint algebra at a quantum level is respecting the symmetries of the theory under 
quantization. It should be noticed that imposing the correct commutation relations between constraints at a quantum 
level without anomalies is far from a mere technicality. For instance in the case of string theory it is what determines 
that the theory of a bosonic string only is consistent in 26 dimensions. 

The canonical formulation of general relativity has constraints. The issue of the commutation relations of the 
constraints therefore appears for any quantized version of the theory that one may consider. In this paper we will 
concentrate on the attempt to quantize General Relativity based on the Loop Representation JIJ constructed from 
Ashtekar's new variables |Q. The issue that we want to address is the consistency of the constraints that have been 
proposed in . Although solutions to the proposed quantum version of the constraints have been found [||| , this 
does not imply consistency of the constraint algebra. 

We will perform the calculation directly in the Loop Representation. This will serve two purposes: firstly it shows 
that the expressions presented in j3|JJ] are well defined in the sense that commutators of the constraints can be 
computed. Secondly it provides an excellent example of the use of several techniques to perform computations in 
loop space. We have made this paper deliberately detailed and explicit, at the risk of being somewhat long, with 
the purpose of presenting in a clear fashion the techniques used. These techniques can be applied to other problems 
concerning the loop representation. 
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The problem of the constraint algebra in quantum gravity has received a fair share of attention, both in terms of 
traditional geometrodynamic variables and in terms of Ashtekar new variables. In particular it is quite clear that the 
issue cannot be completely analyzed on formal grounds: a suitable regularization is needed to avoid ambiguities that 
can arise in the constraints. A thorough discussion of these and other issues, as well as a historical introduction to 
the subject can be found in the papers by Tsamis and Woodard and Friedman and Jack M. 

In spite of what we have just said, in this paper we will only discuss the issue of the formal unregulated commutation 
relations. The motivation is twofold. First of all, not even this calculation has ever been performed in terms of loops. 
Second, our main purpose is more to gain insight on how to operate with the recently introduced expressions for the 
constraints in the loop representation rather than to settle the issue of their consistency. The constraints in the loop 
representation have been used to find solutions to the equations Q, but operating twice with a constraint requires 
considerably more care than originally expected. Several lessons about how to operate with the diffeomorphism and 
Hamiltonian constraints in terms of loops will be learnt during the discussions of the different commutators. A sharper 
understanding of how these operators act is gained through this analysis. 

This paper is organized as follows. In section II we will briefly recall the definition of the constraints of canonical 
general relativity and the constraints of canonical quantum gravity in the Loop Representation In section III we 
compute the commutator of two diffeomorphism constraints, in section IV the commutator of a diffeomorphism with 
a Hamiltonian and in section V the commutator of two Hamiltonians. Several technical aspects of the computations 
are discussed in the appendices. 



II. CONSTRAINTS IN THE LOOP REPRESENTATION 



There exist well established canonical formulations of general relativity. The more traditional one is based on the 
use of a three dimensional metric q a b and its conjugate momentum density Tr ab . Recently, a new formulation due 
to Ashtekar has exhibited interesting properties. It is based on the use of an SU(2) connection A l a as fundamental 
variable and its conjugate momentum density is a set of densitized triads Ef. In both formulations the variables are 
not free but are subject to constraints. These constraints group naturally into a vector and a scalar. In terms of 
Ashtekar's new variables, 



C(N) = J d 3 xN a (x)E^(x)F^ b (x) (1) 
H(M) = f d 3 xM(x)Et(x)E b 3 (x)Fl b (x)e^ k (2) 



where F* b is the curvature of the Ashtekar connection. It is important to remark that the Gauss constraint is not 
mentioned because in the Loop representation it is satisfied automatically for all loop functionals. We refer the reader 
to H for notation and conventions. 

In the loop representation of quantum gravity wavefunctions ^(7) are functions of loops 7. The word loop has 
in this context a very precise meaning: it is an equivalence class of unparametrized curves. The equivalence relation 
is that two curves are equivalent if they give rise to the same holonomy for all connections. This means that two loops 
are equivalent if they only differ by "tails" that immediately retrace themselves. 

Because the theory started with an SU(2) connection (the Ashtekar connection), wavefunctions inherit certain 
SU(2) identities (Mandclstam identities), 

*( 7 ) = f^- 1 ) (3) 
*(7i 72) = *(72 7i) (4) 
*(7i 72 73) + *(7i 72 = *(72 7i 73) + *(72 7i 13 1 )- (5) 

In this representation, the diffeomorphism and Hamiltonian constraints can be written by making use of the loop 
derivative. The loop derivative is the derivative that arises in loop space when two loops that differ by an infinitesimal 
element of area are considered close. The definition of the loop derivative of a function of a loop is, 

zwswr) = «m mi±t±4pi^m (6) 

57— »■ a 

where ^7 is a loop of infinitesimal area a ab basepointed at x and connected to the loop 7 , through a path ttq . We 
will usually consider for practical computations infinitesimal loops formed by a parallelogram of infinitesimal sides u 
and v. 
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Notice that the loop derivative is a path dependent object. Loop derivatives act on any loop or path dependence. In 
particular if one computes the commutator of two loop derivatives there will be a nonvanishing action of the second 
loop derivative on the path dependence of the first loop derivative, leading to a commutation relation of the form, 

[A ab «), A cd (x y )] = A ab «)[A cd ( x *)] = -A crf ( X *)[A ob «)] (7) 

where the brackets denote that the second loop derivative only acts on the loop dependence of whatever is in the 
brackets. 

Another relation of importance is to notice that the loop derivative satisfies a Bianchi identity. In order to define 
it, we need to introduce another derivative in loop space, the covariant derivative. The covariant derivative D a acts 
on functions of paths by appending an infinitesimal element at the end of the path along the coordinate direction a, 

(l + eu a D a )^(^) = ^(7T^ u ). (8) 

With this definition, the Bianchi identity reads, 

Z? a A bc «) + D b A ca {<) + £> c A ab «) = 0. (9) 

The constraints of quantum gravity in the loop representation can be obtained using the loop transform. They 
have been discussed in refs HQ so we only discuss them briefly here to fix notation. 
The diffcomorphism constraint is given by, 

C(N)*(y) = J d 3 xN a (x) j> dy b 5 3 (x - y)A ab (^{ 1 ) (10) 

where N a is a vector field on the three manifold along which the diffeomorphism is taken. This operator has been 
known for quite some time to generate deformations of the loop in the argument of the wavefunction corresponding to 
the diffeomorphism performed. As an intuitive picture of this expression one should remember that in the Ashtekar 
formulation of canonical gravity the diffeomorphism constraint is given by J d 3 xN a (x)Tr[E a (x)F a i,(x)]. The expres- 
sion in the loop representation can be heuristically thought of as a replacement E a — > <f dy a and F ao — > A ao . In this 

paper we will only consider diffeomorphisms that leave unchanged the basepoint of the loops we are using, that is, N 
vanishes at the basepoint. 

The Hamiltonian constraint is, 

iT(M)*( 7 ) = lim| J d 3 xM(x) ^dy"^dz 6 * 3 (x-»)/ e (y, 1 ,0[O»v'A o6 (7*) + O VOI ^A O 6(7*)]*(7). (11) 

In this paper we will omit the factor |. This expression requires some discussion. First of all notice that we 
have introduced a regulator f e (y,y') such that lim e _ >o y') — * S 3 (y,y'). The need for this regulator can be 
seen directly from the translation of the expression of the Hamiltonian constraint in the Ashtekar formulation 
/ d 3 xM(x) J d 3 yf e (x,y)Tr[E a (y)E b (x)F a b(x)] due to the fact that the constraint is quadratic in momenta. The 
expression is naively zero if the loop is smooth (no kinks or intersections) due to the fact that A ob is antisymmetric. 
One can check that when the regulator is taken into account the expression also vanishes provided it acts on diffco- 
morphism invariant functions. Since wavefunctions in the loop representation take values on all possible loops we will 
generically consider a loop with a multiple intersection at a given point (the case with more than one intersection is 
trivially included since the action of the Hamiltonian is local). 

A second point to notice is the appearance of the rerouting operator, O yy i. Acting on a function of a loop, this 
operator has the following effect: it takes the loop formed by taking the original loop between y and y' no matter the 
original orientation and then it adds the portion between y and y' passing through the origin o, no matter the original 
orientation. Of course, the points y and y' must lie on the loop, as they do in the expression of the Hamiltonian. 
The reader may be puzzled at this point: if one cuts and pastes a portion of a loop in the opposite orientation, one 
does not have in general a loop any more. This is true. To complete the procedure for arbitrary y and y' one must 
re-close the resulting loops (see appendix A). However, in the limit where regulators are removed, the points y and y' 
coincide with the intersection points and therefore the action of the rerouting simply amounts to a rerouting of certain 
"petals" of the loop. The "re-tracings" introduced to close the loop are of higher order in powers of the regulator in 
the limit where the regulators are removed. 

The reader can realize that although applying these expressions to a function of a loop seems reasonably straightfor- 
ward, applying them for a second time can introduce complications. For instance: what is the action of two successive 
reroutings? What is the action of the loop derivative on the rerouting? These are the kinds of questions we would 
like to address in this paper. 
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III. DIFFEOMORPHISM ALGEBRA 



We now proceed to compute the commutator of two diffcomorphism constraints. The reader may feel this an 
unneeded computation. After all, if the constraint is the infinitesimal generator of diffcomorphisms, it should, by 
definition, satisfy the correct algebra. Although this is a valid viewpoint, we feel that an explicit calculation is in 
order. Furthermore, the calculation can be viewed as a confirmation that the proposed constraint is equivalent to the 
generator of diffeomorphisms in loop space. The calculation is quite instructive in the sense that it makes a nontrivial 
use of several properties of loop space. This calculation is already present in the literature |j) but we feel a more 
careful and detailed derivation is in order. We find it convenient to explain the techniques used in some detail with 
the example of this calculation in order to be able to be more succinct in the following ones. 

We will start with an auxiliary calculation that will prove useful in the sequel. We will evaluate the action of the 
loop derivative on a diffcomorphism constraint. To this aim we construct the following expression, which holds due 
to the very definition of loop derivative, 

(1 + o^AcdW)) I dy a S(y - x)A ab { 1 ^{l) := 
J j 

dy a 5(y - x)A ab ((S lz o y) v )V(5j z ° 7)- (12) 

In this expression, Sj z is the infinitesimal loop added to 7 and connected trough a path from the origin up to 
the point z. That is, we evaluate the action of an infinitesimal deformation of area a cd acting on the infinitesimal 
generator of diffeomorphisms. We now evaluate the right member of this expression and this will enable us to read 
off the action of the loop derivative on the infinitesimal generator of diffeomorphisms. 

With this aim we expand the right member of ( p^ ) , partitioning the domain of integration, and using the definition 
of the loop derivative to expand \&(#7 o 7). We get, 

cd t 



dy a A ab ( 7 y)S(y - x)(l + <r cd A cd ( 7 *))* (7) + 

+{u a 5(z - x)A ab (Y ) + v a 5(z + u- x)A ab ( 1 z + u ) - u a S(z + u + v- x)A ab ( 1 z +u+v ) 
-v a S(z + u- x)A ab (Y +u+v+a )} (1 + <7 cd A cd ( 7 *))#( 7 ) + 

dy a S(y - x)A ab {{8 lz o 7 )V)(1 + ^^(7*)* (7) (13) 

where the terms in curly braces correspond to the integration along the infinitesimal loop £7 which we take as being 
given by a parallelogram of infinitesimal sides given by u a and v b . 
The last term in this expression can be rewritten as 

dy a S(z - y)0(z, y)[l + a cd A cd ( 7o *))[A ab ( 7 ^)](l + ^A e/ ( 7o z ))*( 7 ) (14) 

where <d(z,y) is a Heaviside function that orders the points along the loop. We will be able to combine the zeroth 
order contribution (in terms of the infinitesimal loop) of this term with the first term in ( [l3|) . It should be noticed that 
the first loop derivative does not act on everything to its right but only on the path inside the second loop derivative 
7q, a fact we denoted by enclosing it in brackets, and that was discussed in detail in the previous section. We now 
consider the expansion of the terms containing infinitesimally shifted loop. They can be expressed with the use of the 
Mandelstam derivative, 

A ab ( 7 r u )=(l + ^c)A ab (7 z ) (15) 
A ab (Y- +u+v )^ (1 + v d D d )(l + u c D c )A ab (Y ) (16) 
5{z + u-x)= (l + u a d a )6{z- x) (17) 

We now expand again ( |l~3| ) 

dy*6(v - x)A ab ( 1 y)(l + ^ d A cd {ll))^{l) + 
+{u a 5(z - x)A ab { 1 z a ) + v a (l + u e d e )6(z - x)(l + u d D d )A ab (^) - 
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u a (l + v c D c )S(z -x)(l + v d D d )A ab { 1 z ) - 
v a S(z - x)A ab (^)}(l + cx cd A cd ( 7 *))*( 7 ) + 

^%-x)e(z j2 ;K d A crf ( 7o z )[A Qb ( 7 ^)](l + a e /A e/ ( 7o 2 ))vI/( 7 ). (18) 



'7 

Of the terms in braces, it can be readily seen that only contributions proportional to u a v b are present, neglecting 
terms of higher order. The other terms combine to give the original expression. We can finally read off the contribution 
of the loop derivative, 



A cd ( 7 *) f dy a 6(x - y)A o6 (7»)*( 7 ) = 
2[d [c S(z - x)A d]b { 1 z ) + 8{z - x)D [c A d]b ( 1 z )]^>( 1 ) + 
dy a e(z, y)S(y - x)A«,(7*)[A„ 6 (7*)]*(7) + 



7 

dy a 8(y - x)A ab ( 1 «)A cd ( 1 z )^(l) (19) 

'7 

With this calculation in hand, it is straightforward to compute the successive action of two diffeomorphisms, 

C(N) C{M) $(7) = (20) 
d 3 wN d {w) I dz c S{w - z)A cd { 1 z ) I d 3 xM b {x) I dy a S(y - a;)A a6 (7»)tf (7). 

J -y J J ~y 

Expanding this expression, we get six terms, 

d 3 w J d 3 xN d (w)M b (x) j dz c S{w - z)d c S(z - x) A db (7^)^(7) - 



d 3 w J d 3 xN d (w)M\x) j dz c S{w - z)d d 8{z - x)A cb ( 1 z )^{ 1 ) + 
d 3 w J d 3 xN d (w)M b {x) j> dz c S{w - z)5(z - x)D c A db { 1 z )^{ 1 ) - 
d 3 w ( d 3 xN d (w)M b (x) <f dz c 5(w ~ z)5{z - x)D d A cb (Yo) 1 S>h) + 



j j 7 
d 3 w [ d 3 xN d (w)M b (x) <f dz c <f dy a S(w - z)G(z,y)6(y - x)A cd ( 7o z )[A afc (7^]*( 7 ) + 

d 3 w f d 3 xN d (w)M b (x) I dz c I dy a S(w - z)S(y - x) A ab ( 7o «)A cd ( 7o 2 )*( 7 ) (21) 

We should now subtract the same terms with the replacement N <-» M. Since the calculation is tedious but 
straightforward we describe in words how the terms combine. The fifth and sixth terms, when combined with the 
similar terms coming from the substitution N <-> M cancel taking into account the commutation relations for the loop 
derivatives (0). The first and third term, combined with the first of the substitution N «-+ M form a total derivative. 
The fourth term, combined with the third and fourth of the substitution TV <-> M cancel due to the Bianchi identities 
of the loop derivatives. Finally, the second terms combine to produce exactly C(C^M), which is the correct result of 
the calculation. 



IV. COMMUTATOR OF A DIFFEOMORPHISM WITH A HAMILTONIAN 



This calculation will teach us about two important new ingredients that were not present in the previous calculation: 
the action of reroutings in the Hamiltonian and how to deal with regularization. There has always been concerns about 
the use of a background dependent regularization in the Hamiltonian constraint, since it had potential to interfere with 
the action of diffeomorphisms. With the help of this calculation it is possible to detect the regularization problems 
explicitly, by showing where it is needed to remove the regulators for the computation to work. The commutator is, 
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[C(M)H(N) — H(N)C(M)] = J d 3 wM h (w) J d 3 xN(x) x 

<f dz a 8(z - w)A ab (Y ) <f dyl c <f dy'^5(x - y')f £ (y - y')[O vy , + O yoyl ]A cd (l v )^{l) ~ 

J L J 7 J 7 

/ dy^ <f dy'^S(x - y')My - y')[O vy , + yoy ,}A cd (^ o ) <f dz a 5(z - w)A ab { 1 ^( 1 ). 



(22) 



We will now proceed as in the last section; first, we evaluate the action of the loop derivative of the diffeomorphism 
constraint on the loop derivative and the rerouting operators corresponding to the Hamiltonian constraint. We will 
prove that these terms minus the analogous ones coming from the action of the Hamiltonian on the diffeomorphism 
cancel each other. The remaining terms, which are the consequence of the action of the loop derivative on the loop 
dependence of integrals will give rise to the correct commutator. Since the terms involving O yoy > behave in an 
analogous fashion to the ones that depend on O yy i, we will only concentrate on these. 

The first contribution of the action of the diffeomorphism on the Hamiltonian is, 



<f dz a 8(z -w)<f dy [c <f dy ,d ^(x - y')f e (y - y') 

J 7 J 7 J 7 



A ab (Y )[O yy ,A cd ( 7 «)*(<y)]. 

We now separate the outer integral into two portions, using the same notation as in the rerouting operators, ie. 
lyoy' would be the portion of the loop 7 that goes from y to y' through the origin. We also perform a split of one of 
the inner integrals, 



<f dy [c I dy' d ^ I dz a + I dz a 
+ l dy( c I dy' d ^ I dz a + I dz a 

J 1 J'Yvo lJ-y„„i J~t„' n „ 



S(z -w)5(x -y')f e (y -y') x 

*A ab { 1 z )[O vyl A cd { 1 y)*{ 1 )] + 
5(z - w)S(x - y')f e (y - y') x 

xA ab (Y )[O yy ,A cd (<yy o Ml)}- 



(23) 



After commuting the rerouting and the derivative (see appendix A), and taking into account the sign introduced 
in this process we get, 



dy 



id] 



lay 



ly', 



+ l dy^ c I 



dy 



id] 



dz a + 



dz a 



dz a 



ly o y' 



dz a 



TV, 



8(z - w)5(x - y')f € (y - y') x 

x0 !/!/ ,A a6 (7 o 2 )[A cd (7^*(7)] 
6(z - w)5(x - y')f e (y - y') x 

xO, J/ ,A a6 ( 7o z )[A cd (7^)*(7)]. 



(24) 



The following step consists in noticing that the minus sign in front of the integral can be absorbed by integrating 
along the loop in the opposite direction, which we denote with an over-bar, 



+ & dy* 



dy 



id] 



lay 



dz a + 



dz a 



ly o y' 



6{z- w)S(x - y')f e (y - y') x 

xO^A a6 ( 7 *)[A C cK7 y M7)] + 



dy 



id] 



lyc 



dz a + 



dz a 



5(z - w)5(x - y')f e (y - y') x 

xO^A a6 ( 7 *)[A cd (7 y M7)], 



(25) 



which in turns allows us to recombine the terms again into a single loop integral, making use of the definition of the 
rerouting operator, 



G 



di c f dy'^Oyy, f dz a S(z - w)S(x - y')f e (y - y') x A afc ( 7o z )[A cd ( 7 ^)*( 7 )] + 
+ <f di c <f dy^O yy , <f dz a 5(z - w)8(x - y')f e (y - y') x A afc ( 7o 2 )[A cd ( 7 ^)vI/( 7 )]. (26) 
So the final result of this manipulation can be written as, 

dy^ c <f dy'^Oyy, <f dz a S(z - w)S(x - y')f e (y - y') x A afc ( 7o 2 )[A cd ( 7o ^( 7 )] = 
dyl c (f dy'^5(x - y')My - y>) x O y y , 

7 Jl 

dz a S(z - w)[A cd (^ )A ab (Y ) + e(z,y)A afc ( 7o z )[A crf ( 7o ^)]]*( 7 ). (27) 

'7 

The corresponding terms arising from the action of the Hamiltonian constraint on the diffeomorphism are, 
dy^ I dy'^Oyy, I dz a S(z - w)5(x - y')f e (y - y') x A cd ( 7 ^)[A afc ( 7o 2 )vI/( 7 )] = 
= <f dyl c <f dy'^Oyy, <f dz a S(z - w)6(x - y')f e (y - y 1 ) x 

J -y J -y J -y 

x[A afc ( 7o z )A cd ( 7 «) + e(y^)A crf ( 7 ^)[A afc ( 7 2 )]]vI/( 7 ). (28) 



When we subtract between (|28|) from ( |27| ) one immediately notices that the net result is zero if one uses the 
commutation relation for the loop derivatives (Q). 

We now concentrate on the terms that will give rise to the result of the commutator. These terms appear when 
the loop derivative acts on the loop dependence of the integrals, exactly as in the commutator we considered in the 
previous section. 

The contribution from C(M)iJ(7V)*( 7 ) is 



d 6 wM b (w) j dz a 6(z - w) J d 4 xN{x) x 

x{2^ dy' d 8{x - y')Df a [f £ (z - y')O zy , A^l)]^^) + 

+2 j> dy c Df b [S(x - z)f £ (y - z)O y ,A o]c ( 7 *)]*( 7 )) . (29) 
And the contribution from H(N)C(M)^(-/) is 

d 3 xN(x) f d 3 wM b {w) <f dy^ c <f dy' d ^5{x - y')f e (y - y') x 

J J -y J <-y 

x Uo yy .Dl[5(y - w)A d]b ( 7 y)M 7 )\ . (30) 

The result of computing the commutator is, 

2 I dz a 8(z -w)<f dy' d 8(x - y')Df a [f £ (z - y')O zy , A b]d ( 7o z )]vI/( 7 ) + (31) 



+2 f dz a 6(z -w)f dy c Df b (S(x - z)f e (y - z)O y 2 )A a]c ( 7 «)v&( 7 ) 

J 'y J -y 

j 7 



dy^ c f dy'^S(x - y')Uy - y')20 yy , D\ c [5{y - w) A d]b ( 7 ^)]M/( 7 ). 
We now explicitly write the antisymmetrizations of the first two terms, and use Leibnitz' rule in the last one we get, 
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d A wM b {w) J d A xN{x) x 
x(<f dz a S(z -w)<f dy' d 6(x - y')D*[f £ (z - y')O zyl A bd {l z Q W{l) ~ 

dz a S(z -w)<f dy' d S(x - y')D z b [f e (z - y')O zy ,A ad (j*)}*(j) + 
+ f dz a S(z-w)(f dy c D z b {8{x-z)f e {y~z)O yz )A ac {lt)^{l)- (32) 



dz a 5{z -w)f dy c D z a (S(x - z)f e (y - z)O y ,)A bc ( 7 £)*( 7 ) 

7 J~f 



dyl° f dy'Mx - y')My - y')O y - w)A db (^ ) - d y d 5{y - ™)A cb ( 7 ^)]4>( 7 ) 

dy^ c I dy'^8{x - y')My - y')O yy> 8{y - w)[D v c A db {~f v ) - £>»A c6 ( 7 »)]* ( 7 )' 



'7 J l 

The first and fourth terms can be integrated by parts on z. Taking into account that the diffeomorphisms we are 
considering have trivial action at the basepoint of the loops, the resulting terms cancel with the fifth term. In order 

for this cancelation to occur we have to remove the regulator so that d v d 8{y' — w)f e (y — y') — d%6(y — w)f e (y — y'). 
The remaining terms are the second and third. By using the Bianchi identity in the last term of (p3), this equation 
may be rewritten as 



d 3 wd 3 xM b (w)N(x) x 

dy c S(y -w)j> dy' d 5{x - y')Dl\f e {y - y'^A^)}^) + 

+ <t dy c S(y -w)<f dy' d Dl(6(x - y)f e (y' - y)O y , y )A cd (l V o'Ml) + 

+ j> dy^ j> dy' d Xv ~ w)5(x - y')f £ (y - y')O y ^A od ( 7 *)*( 7 )) . (33) 
By using Leibnitz rule in the first term we get, 
J d 3 wd 3 xM b (w)N(x) x 

x ( - j> dy c S(y -w)j> dy' d S(x - y')Dl[f t {y - y')O yy> ] A cd ( 7o ^)*( 7 ) + 

+ j dy c S(y -w)j> dy' d Dt(S(x - y)f e (y' - y)O y , y )Kd(rf)m)) ■ (34) 



'7 •>! 

Taking into account the definition of the rerouting operator, we have 



D y b lMy - y')O yy> ] = Dlf e {y - y')O yy , = -D\ f e (y - y')O yy . (35) 



and using the Leibnitz rule again in (34) we get 

d 3 wd 3 xM b (w)N(x) x 



x(+<fdy c <f dy' d Dt[5{y - w)6(x - y')f £ (y - y')O yy> A cd ( 7 *)]M>( 7 ) + 



+ 



I dy c I dy' d Dl[5{y - w)5(x - y)f € (y' - y)O y , „A cd ( 7 »')]*( 7 ) + 



dy c S(y -w)i dy' d d y b S(x - y')f e (y - y')O yy , A cd {l v )*{l) + 



7 J l 



7 J l 



dy c d v b 6(y -w)(f> dy' d 5(x - y)f e (y' - y)O y , 9 A ci ,( 7 «')*( 7 ). (36) 
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The second term cancels with the first after interchanging y and y' (and again removing the regulator so that 
S(y — w)f e (y — y') = S(y' — w)f t (y — y')). If we interchange y and y' in the last term and integrate by parts in x and 
w we get the final result, 

f d 3 wd 3 x(N(x)d b M b (w) - M b (w)d b N(x)) I dz a I dy' d x 

J J -y J -y 

xS(z - w)S(x - y')f e (z - y')A ad (Y )O zy ^h) = H(C^N) (37) 
which reproduces the classical commutator. 

V. COMMUTATOR OF TWO HAMILTONIANS 

The commutator of two Hamiltonians introduces another new computational requirement. While computing the 
commutator, one gets the product of two rerouting operators: one per each Hamiltonian. The result one expects 
involves only one rerouting, the one due to the presence of the metric in the function smearing the resulting diffco- 
morphism. We will give greater details of this aspect of the calculation in appendix B. 

The commutator of two Hamiltonians can be written as, 



[H(N)H(M) - H(M)H(N)]V(L) = J d 3 xN(x) J d d wM(w) x 
/ dy^ c <f dy'^6(x - y')U(y - y')[O vy , + CW]A cd ( 7 «) x 
dz^ I dz"^6(w - z')f p (z - z')[O zz , + O zo ^]A pft ( 7o z )*( 7 ) - 
J d 3 xN(x) J d 3 (w)M(w) x 

dz^ <f dz" l ^S(w - z')f p (z - z')[O zz > + CW]A p/l ( 7o z ) x 
I dy^ c <f dy' d ^6(x - y')f £ (y - y')[O yy> + yoy ,]A cd (l v o )^{l)- (38) 

We now proceed as in the previous case. Since several terms in the commutator cancel in similar fashion, we will 
show the cancellation explicitly for one case and will indicate how to proceed with the others verbally. 

Let us start by considering the terms that arise from the action of the derivative operator of the first Hamiltonian 
on the integrand of the second operator. We will discuss later on the terms that arise from the action of the derivative 
operator on the loop dependence of the integrals of the second Hamiltonian, which will give rise to the nonvanishing 
part of the commutator, which reproduces the classical result. 

The calculation will be performed as follows. We will consider each of the above mentioned terms in H(N)H(M) and 
perform operations that show that they are actually term in H{N)H{M) and therefore cancel in the commutator. 
This will involve proving identities among integrals along petals of intersecting loops. For this it is convenient to 
introduce loop diagrams. In these diagrams we represent each loop as an oriented segment (from left to right) and 
we mark the ordering of the different variables in the integrands along the loops. We will see that in several cases 
the order in which variables appear implies that the integrals are actually zero. Several operations we will perform 
involving the rerouting operator will resemble the ones we performed in the previous section. 

We now consider one of the contributions indicated above, extract the rerouting operator outside the loop integral 
as we did in the previous section and we get, 



(fdj/ c (f dy' d ^(x-y')f e (y-y')x 

I dz^i dz' h ^{w-z')f p {z-z')O yy ^Ul V o)[O zz '\ h {n^{l)] 

J[o yy n] oy > Ao yy n]yo 

+ <f dyl c <f dy"^S(x - y')My - y') x 

J-y Jly o 
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d7> <f dz' h ^5{w - z')f p (z - z')0 y0y ,A cd ( 7 y)[0 zz ,A ph (y z )*(j)} = 

[o yoy n] oy > J[o yoy n] V o 



= fdy^f dy' d ^(x-y r )f e (y~y')x 

J J J toy 

dz^ <f dz^S(w - z')f p (z - Z^Oy^Ozz^M^M^il)] 

\o yy n] oy > J [O yy n]yo 



-fdy^f dy' d ^S(x - y')My - y') x 

dz^ I dz' h ^S(w - z')f p (z - z')O y o,yO Z2 'A C(i (7 o y)[A^( 7o z )^(7)]. 

[o yoy n] oy > J[o yoy n]yo 

We now use the commutation relation for loop derivatives (M) to obtain the identity, 



(39) 



A cd (7^)[A p ,( 7o z )*(7)] = 

= [A p ,( 7o z ) A crf ( 7 *0 + 9(z - y)A C(i (7^)[A p ,( 7o z )]]*(7) - 

= [A p ,( 7o 2 )[A C(i (7 y )] + A cd ( 7 ^)A p/l (7 2 ) - e(y,z)A p ,(7 z )[A cd ( 7 ^)]]*(7) = 

= [e(z., y )A ph (Y )[A cd ( 1 y)} + A cd ( 7 »)A ph ( 7o a )]*(7) = 

= A p/l ( 7o 2 )[A C£i (7^)*(7)]. (40) 



It is now useful to draw the loop diagrams corresponding to the integrals present in (|39j). In the diagrams we 
indicate the relative domain of integration of the different variables in the first term of the expression, 

1 — I — I — I — I 1 — I — I — I — I 

z y' y z ' o z ' y y' z 

By looking at the diagram, we see we can rewrite the integrals in the following way 

1 dz^ <f dz'^5{w - z')f p (z - z') x 
dy^ c I dy'^S(x - y')f e (y - y')O yy ,O z z ,A p/l (7^,(7^(7)] - 



w — z') f p (z — z') 

J T J 7„ > 



dy^ c I dy^5(x-y')My-y')O yoy/ Ozz'A ph (Y )[A cd (l v )yh)]- (41) 

We now use the fact that O yy >0 Z z' = Z z'O yy ' and commute rerouting and derivative operators again. 
& dz [ P & dz ,h] S(w - z')f p (z — z') X 

dy^ c I dy'^S(x - y')f e (y - y')0,,'0,^ A^( 7o 2 )[A cd (7 y )f (7)] - 



<b dz^ <h dz' h ^8{ w — z') f p (z — z') 
J 7 Jjoz 



dy^ j dy^5(x - y')f e (y - y')O z z ,O yoy ,A ph ( 7o *)[A cd ( 7 *0* (7)] 

= <f> dz^ i dz"^5{ w - z')f p (z — z') X 

/ dyl c I dy'^S(x - y')f e (y - y')0 Z7 j A p , ( 7o 2 ) [O, y A cd ( 7o y )*(7)] 4 
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+ <b dzb <p dz" l ^S(w - z')f p (z - z') x 



lo 



dyl c f dy' d ^(x - y')My - y')O z z , A^Dpy oy , A cd {l v )*{l)] 

[o z ,n],*> J[a*,H]*v 



w — z')f p (z — z') 



J 7 J 7 

dyl c / dy'^S(x - y')U{y - y')O zzl A ph (Y )[[O yy , + O yoy ,]A cd {ll)^{l)]- (42) 

This proves that all terms with four loop integrals that appear in the commutator ( |38| ) cancel among themselves 
(each order of the operators in the commutator produces an equal contribution). 

We now need to address the terms with three loop integrals. In order to do this we begin with the following 
observation. If in the second term of expression ( |38|) we rename pairwise the names of the following variables, z and y, 
z' and y', p and c and finally h and d. The only real change in the last term of the commutator is that x gets replaced 
by w in the delta's arguments. Then, we may work with the first term and only at the end add the contribution of 
the second, which is similar in form. As before, we need only consider the proof for O yy iO zz i. Since in the terms we 
are considering the loop derivative acts only on the loop dependence of the integrals, terms with different reroutings 
can be treated in a similar fashion. 

We begin studying the action of the loop derivative of the first Hamiltonian on the third and fourth loop integrals, 



dzb <f dz'VS(w - z')f p (z - z')O zz ,A ph ( 7 *)*fr) - 
dz^ f dz^6(w ~ z')f p (z - z')0 zz ,A ph (Y o )^h) + 



+ J> dz^ f dz lh ^{w - z')f p (z - z')0 Z2 'A p/i (7 o z )f ( 7 ) 

- I d?> <f dz' h ^(w ~ z')f p (z - z')0 zz ,A ph (Y o )^(l). (43) 

In the first integral we need to do the following substitution (valid only in the limit when the regulator is removed), 

f P (z-z')O zz ,A ph (Y )*(l) = -f p {z-z')O z , oz A ph {^). (44) 

As an aid to understand this expression one should remember that the corresponding expression in the connection 
representation would be, 

f p (z - z')Tr[F ph (z)U zz ,U zoz ,] = -f p {z - z r )Tr[F ph {z')U z , oz U z , z ] (45) 

Finally, working in an analogous way as in the previous section the result is, 

dyl c I dy' d ^(x - y')My - y') x 



7 j 1 



O yy , (tfdz'Hw - z)Df c [f p (y - z)O zoy A d]p (Y )Ml)+ 

+2 / (WDmw - y)f P (y - z)O zy A d]p ( 7 z )]*(j)) • (46) 



Considering the fact that the Mandelstam derivative commutes with the loop integral and the rerouting operator, 
as it can be seen through the results of appendix A, we get, 

dy c (f dy' d 5(x-y')My-y') x 

7 J 7 

Dl\O vy , I dzP[S(tv - z)f p (y - z)O zoy + S(w - y)f p (z - y)O z y ] A dp (Y )}* (j) - 
D y d [O yv> j dzP[S(w - z)f p (y - z)O zoy + S(w - y)f p (z - y)O z y ] A cp ( 7 *)]*( 7 )) , 

(47) 
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and removing the regulator we get, 



dy c j> dy' d S(x-y')My-y') x 

7 Jl 



d«S(w - y)O vy , <f dz*>[f p (y - z)O zoy + f p (z - y)O zy }A dp ( 1 z )}^( 7 ) + 

J y 

+8(w - y)DV[O vv , <f dz*>[f p {y - z)O zoy + f p (z - y)0^]A dp ( 7 *)]*( 7 ) - 

J-y 

-d v d S(w - y)O yy , <f dz"[f p (y - z)O zoy + f p (z - y)O zy }A cp ( 1 z )}<f( 1 ) - 

J y 

~8(w - y)D y d [O vv , j dzP[f p (y - z)O zoy + f p (z - y)O zy ] A cp ( 7 *)]*( 7 )) (48) 

In the limit in which the regulator is removed the product S(x — y) x5(w — y) is symmetric in x and w and therefore 
some of the terms in the above expression cancel with those coming from the other term in the commutator. The 
remaining terms are 

<f dy c I dy' d S(x-y')f e (y-y>) x 

J y J y 

(-8?5(w - y))O yy , I dzP[f p (y - z)O zoy + f p (z - y)O zy }A dp (Y;)}*h) - 

J y 

- <f dy c <f dy' d 5(x-y')f e (y-y') x 

J y J y 

(-d$5(w - y))O yy , I dzP[f p (y - z)O zoy + f p (z - y)O zy ] A cp ( 7o z )]M/( 7 ). (49) 

J y 

The following step is to integrate by parts 
J d 3 xN(x) J d 3 wM(w) x 
x(^jdy c ^dy' d S(x - y')f £ (y - y') x 

(-d™5(w - y)) O yy , <f dzP[f p (y - z)O zoy + f p (z - y)O zy ] A dp ( 7o z )]*( 7 ) - 
- I dy c I dy' d 5(x-y')f e (y-y') x 

J y J y 

{-d$5(w - y)) O yy , j dzP[f p (y - z)O zoy + f p (z - y)O zy }A cp (Y )}*(^ = 
= - J d 3 xN(x) J d 3 w D™ (m(w) j dy c j> dy' d 6(x - y')f £ (y - y')S(w - y) x 



'7 J 7 

Oyy, i d Z Vf P (y - Z) [0 ZO y +0 Z y]A dp { 1 l)]^( 1 ) ) + 



+ f d 3 xN{x) f d 3 w d™M{w) x I dy c I dy' d 5(x - y')f e (y - y')S(w - y) x 

J J J y J y 

Oyy, <f dz"f p (y - Z )[0 ZO y+ O Z y] A^)^ (j) + 

+ f d 3 xN(x) [ d 3 w D d (m(w) I dy c <f dy' d 5{x - y')f e (y - y')S(w - y) x 

J J \ J y J y 

Oyy, l d Z Pf p (y - Z )[O ZO y+ O z y } A cp (<£ ) ] * ( 7 ) ^ 



(50) 
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d 3 wd%M(w) x &dy c & dy' d 5(x - y')f e (y - y')5(w - y) x 



-J d 3 xN(x) J , 
O vy , <f dz"f p (y - z) [O zoy + O zy }A cp ( 1 z )}*( 1 ), 



(51) 



and to notice that since we are considering only the case of compact manifolds, the total derivatives that appear in 
the first and third terms do not contribute. The resulting terms can be rewritten as, 



(52) 



d 3 wN(x) d™M{w) x2<f dy^ I dy' d ^6(x - y')f t (y - y')5(w - y) x 
O yv , j> dzPf p (y -z)[O zoy + zy ]A dp { 1 z o )]^>{ 1 ). 



Taking into account the discussion at the beginning of this section, it is very easy to see (due to the symmetry in 
changing x and w in the deltas) that the other term of the commutator corresponding to the last one obtained is, 



-J d 3 x J d 3 wM(w)d x c N(x)x2^di c £dy' d ^(x-y')f e (y-y')6(w-y) 
O yy , <f dz*f p (y - z) [O zoy + O zy ]A dp ( 1 z )}*( 1 ). 



(53) 



The result then is, 



J d 3 x J d 3 w (N(x) d™M(w) - M(w) d*N(x)) x 
( ldy c l dy' d S(x - y')f £ (y - y')S(w - y) x 

\ J "y J 'v 

Oyy, l d Z P f p (y ~Z)[0 ZOV + Z y] A dp ( 7 * )] * ( 7 ) - 



dy d & dy ,c S(x - y')f € (y - y')6(w - y) x 

Oy y* j dz? f p (y-z)[O ZO y+O z y]A dp ( 7 *)] * ( 7 )^ . 



(54) 



We now make some changes in this last expression so that we can use the equality of appendix B. Noticing that, 



S{x - y')fe(y - y')S(w -y) = 5(x - y')f € (w - x)S(w - y) 
we can rewrite the expression of interest as, 

d 3 x [ d 3 w (N(x) d™M(w) - M(w) d*N(x))f e (w - x) x 



(55) 



and then as, 



J dj 3 5(w -y) f dz 3 f p (w -z)<f dy c 6(y -y) <f dy' d S(y' - x) x 

J J J -y J -y 

Oyy, <f dz"6(Z -Z)[0 ZO y + Z y] A dp (<£ )] * ( 7 ) - 

-O y ,yjd Z P5(Z-Z) [0 ZO y, + Z y,]A dp { 1 Z )]^( 1 ) 



J d 3 x J d 3 w (N(x) d?M(w) - M(w) d*N(x))f e (w - x) x 

J dy 3 S(w - J/) J dz 3 f p {w -t)j> dy c S(y - y) jf dy' d S(x - y 1 ) x 

<f dZ P 5(Z - Z) [Oyy,[0 ZOy + Z y}+ O >y, y[0 Z y, + Z y, ] ] A dp * ( 7 ) . 



(56) 



(57) 
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Finally, using the result of appendix B we can rewrite the above expression in terms that only involve one rerouting 
operator, as is required by the desired final result (which only involves one rerouting coming from the metric that 
appears smearing the diffeomorphism constraint). As a consequence of this the commutator is, 

d 3 x J d 3 w(N(x)d™M(w) - M(w)d*N(x))f £ (w - x)x 

dy 3 6(w -y) dz 3 f p (w - z) x 



(58) 



dy c 6(y -y)i dy' d S(x - y 1 ) f dz p S(z - z) [O^A^l ) - O y , 2 A ip ( 7 J)]*( 7 ) 
+ j> dy c S(y -y)j> dy' d 8{x - y') X j> dz*5(z - z)[^A dp {^) - \^{l v )Mi) 



We now prove that in the above expression two of the terms cancel each other and the other two give rise to the 
desired result. Consider the following two terms, 



dy c S(y -y)i dy' d S(x - y') d> dz?5(z - z) O y , ^(7^(7) 



O 1 t 



+ <f dy c 5(y -y)cf dy' d 5(x - y') <f dz p S(z - z)±A dp ( 7 %)y(y). 



J 1 



(59) 



We can rewrite them as, 



d 3 y'5{x - y>) f dy c S(y -y) i dy' d S(y' - y') f dz p S(z - z) O y > z A dp ( 7 *)*( 7 ) 
+ J d 3 ^S(x -tf)j> dy c S(y -y)j dy' d S(W - y') j dz p 6(z - z^A^*)* (y) . 



(60) 



It is convenient to write the product of the last two loop integrals symbolically and study them with diagrams, 
assuming for example that y' precedes y in the loop 7. The alternative case can be studied in a similar way. So, 



(61) 



/ dy' d 




J dz p + 


f dz p 


'fay 


L 7o y' 


"Yy y' 


Jlyo 



z y 


y y' 


z y z 


f y 




h / dy' d 


J dz p 4 


J dz p + 


L 


dz p 


•>ly° 




~^y y' 


10 y' 




z y y' 


y z 


y' y y' 


z 




= f dz p 


[ dy' d 


+ / dy' d + 


•f~oy 


dy' d 


J toy 


- ->lzy 








+ f dz p 




+ / dy' d H 


•f-Yoy 


dy' d 


Jjya 


'2 y 









(62) 



(63) 



We can then rewrite expression (|60|) as, 

1K? 



d 3 y'6(x - y>) (f> dy c S(y -y) f dz p S(z - z) f dy' d 5{y' - y') O y , z A dp (H)*{y) 

7 J 7 JO Z y1 

f / d 3 ^S(x -J)j> dy c S(y -y) <f dz p S(z - z) <f dy' d S(y' - 7)^(7^(7), 



(64) 
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and making the change y' by z (and removing both regulators so that one can change y' by z only in the last two 
deltas) we get, 



/ 



d 3 y'6(x - y>) f dy c S(y -y) f dy' p S(y' - y') f dz d S(z - z) O zy , A rfp ( 7 ^)*( 7 ) 



+ J d 3 y'5(x - y') j> dy c S(y -y)j> dy lp 5{y' - y') j> dz d S(z - z)iA dp ( 7o ^)*( 7 ). (65) 
Finally, taking into account that O y > z Ad p ( 7 ^) = — O z y > Ad p ( 7 ^) , and changing d by p we get, 



d 3 y'S(x - y>) f dy c S(y -y) f dy' d 6(y' - y 1 ) f dz?6(z - z) O y , z A dp ( 7 ^)*( 7 ) 
d 3 ^6(x -y~i)j dy c S(y - y) j> dy ,d 6(V - y') j> d#5{z - z)^ A dp ( 7o a )*( 7 ). (66) 

We have therefore proved that expression |6^ is equal to minus itself. Using the same tools one can demonstrate 
that the two remaining terms in ^8] are equivalent to, 

d 3 x J d 3 w{N{x)d™M{w) - M(w)d*N(x))f e (w - x)x 

dy 3 S(w - y) I dz 3 f p (w -z)x 



(67) 
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dy c S(y -y)f dy' d 8{x - y') f dz p 8{z - z) O y y A dp ( 7 ^)*( 7 ) 



+ j> dy c 5(y -y)j dy' d 5{x -y')xj d Z PS(z - z)± A dp ( 7o z )vI/( 7 ) 



Inserting the factor | that was omitted, this last expression is the regularized form of the product of the metric 
operator times the diffeomorphism constraint, which is consistent with the classical result (see appendix C). 

VI. CONCLUSIONS 

In this paper we have studied the formal commutation relations of the constraints of quantum gravity in the loop 
representation. We have shown that the expressions available for the constraint are operationally useful to compute 
the commutators. We use a background dependent regulator and only recover the correct commutation relations in 
the formal limit in which the regulator is removed. Our computations set the stage for a future computation keeping 
the regulators at the lowest order. 

VII. ACKNOWLEDGEMENTS 

This work was supported in part by grant NSF-PHY93-96246, NSF-PHY-92-07225 and by research funds from the 
University of Utah and The Pennsylvania State University. We also acknowledge support of Pedeciba and Conicyt 
(Uruguay). 

APPENDIX A: 

We need a gauge invariant regularization prescription. If one performs a point-splitting at intersections, one needs 
to close the loops to preserve gauge invariance after the rerouting process. In the following example we denote the 
intersecting point as I. We will use the following notation in this appendix: ^ yy i denotes the portion of 7 from y 
to y' without passing through the origin no matter the original orientation of that portion. The other cases that 
appear have an analogous meaning. Explicitly, the gauge invariance-preserving action of the rerouting operator can 
be defined as, 
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f e (y - y') O yyl *( 7 ) = f e (y- y')^{O yyn ) = 

fe(y-y') ^{lyy' ° ly' I ° ll y ° ly oy' °ly'I°Hy)- 



(Al) 



Let us now consider the action of the rerouting operator on the loop derivative. Consider that the point z is in the 
piece Ju'oy {y before y', for example) Then, 



O yy . A dp (Y Ml) = -Ad^Oyy,*^) = A dp ([O yy ,jr o )*(O yyn ). 



(A2) 



The minus sign appears because the rerouting affects the loop dependence of the loop derivative. It is a situation 
analogous to action of the parity operator in quantum mechanics on the partial derivative, 



Pd x V{x) = -d x PV(x) = d x ^(x') 



(A3) 



If the point z of the loop derivative is in "f yy i there is no sign change. 

Let us now consider specifically the case that appears in the Hamiltonian constraint, where the loop derivative is 
evaluated at one of the points of the action of the rerouting operator, 



v dp fe(y - y') Oyy> A d p( 7 ^)*(7) = f e (y - y') 5y^(j yy ' °7s'J 7/v°7W °ly'i °Hy) 
where the result of the operator 8 y 

S y^(lyy' ° ly' I ° H v ° ly oy' °ly'I°Uy) = ^{h°lyy' ° ly' I ° H y ° ly oy' °ly'I°Uy) 

y, O Jy, J O 7/ y O "f yoy l O Jy/ J O Jjy) 



(A4) 



(A5) 



Notice that the definition of the operator 8 y implies that the action of the loop derivative is with respect to 7 and 
the rerouting operator acts afterwards. 



APPENDIX B: 



In this appendix we will prove that expressions (|57]) and (|58|), ar e equal jLOj. First we write the Mandelstam 
identities for wavefunctions. For loops 7,71,72, and 73, 



$( 7 ) = *( 7 -i) 

*(7i 72) = ^(72 7i) 

*( 7 1 72 73) + *( 7 1 72 l^ 1 ) = *( 7 2 ll ° 73) + *( 7 2 7l 73~ )• 



(Bl) 
(B2) 
(B3) 



Considering equation (A5), (with p = z) we can write the action of the loop operators in expression (|5 

S Z ^(Oyy'0 ZOy -f) +S Z ^{Oyy'0 Z yj) +S Z ^(Oy'yO ZO y'-f) +S Z ^(Oy'yO Z y'-/) 



(B4) 



Additionally, using the results of the appendixes A and the Mandelstam identities of Appendix B we can show 
(p is the intersecting point of 7 i 72 73), 



<5 P *(7i 72 73) + <5 P *(73 7i 72) + <^ P *(72 li ° 73) + <S P *(7i ° 73 ° 72) = 
2[<5 P *(73 72 7i) + ^<V*(73 °72 7i) ~ s P ^(li 73 72) - ^<V#(7i 73 °72)]- 
Using this last identity, the integrand in ( |57| ) can be written as, 

2[S y ^(O zyl ) + ~<V*(7) - S y *(O y , z7 ) - i* v *( 7 )]. 



(B5) 



(B6) 

On the other hand, we can make a partition in the loop integrals so that expression (^7|) splits in 6 terms. Symbol- 
ically we can write, 



dy c f dy' d f dz p 

■'7 '0„'1 



dy c 



dy 



id 



dz p 



dz p + 
Ol ov ' Jo lv 



dz p 



dy 



id 



lyo 



(B7) 
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Then, combining (B6) and (B7), we finally get expression (pq), which in shorthand writing is, 

2(fdy c (f dy' d <f dz p [O zy & dp {^) - O y , z A dp ( 7o a )] + 

J-y J-y JO yy n 

+2 j dy c j dy' d j> dz^ dp {l y o) ~ (B8) 

As an aid to visualize the meaning of the identity (B5), its counterpart in the connection representation is the 
following identity between SU(2) matrices, 

Tr[x d BA^C] + Tr^CBA- 1 } + Tr[x d A~ l B~ l C] + Tr{x d B~ 1 CA~ 1 ] = 

2[Tr[Bx d C] Tr[A] ~ ^Tr[ABx d C] - Tr[Ax D B] Tr[C] + ^Tr[Ax D BC}} 

(= 2[^Tr[ABx d C} + Tr[Bx d CA~ 1 } - ^Tr[Ax D BC] - Tr[Ax D BC- 1 ]]). (B9) 

Where x d = \a d , and a d 1 < d < 3 are the Pauli matrices and A, B and C, are complex SU(2) matrices. This 
identity can be easily proved noting that, 

Tr[x d PQ~ 1 } + Tr{x d Q- 1 P] = 2 Tr[x d P] Tr[Q] - Tr[x d PQ] - Tr[x d QP] (BIO) 

and, 

Tr[p- 1 x d Q] = -Tr[x d P] Tr[Q] + Tr[x d PQ] (Bll) 

Tr[Qx d p- 1 ] = Tr[x d QP] - Tr[x d P] Tr[Q]. (B12) 



APPENDIX C: 

To obtain the expression for the metric operator in the loop representation we first compute it in the connection 
representation and apply the loop version. The action of the metric operator is, 



f e (w-z) — —rv - .... Tr[Ufr)] = f € (w- z) & dz a 8{z-z) 



8A{w) k b 6A{z)l 



7 



dw b 5{w-w)Tr[U{ lwz )U^ woz } + f dw b S(w~w)Tr[U( lzw )U(j zo , 



+Mw - z)~ j dz a S(z dw b S{w - w)Tr[U(j)} 



dw b S(w -w)O wz Tr[U(j)} 



z ~fe{w ~ z) f dz a 5(z - z) 
-f e (w - z)i j> dz a 6(z -z)<j> dw b S(w - w)Tr[U(-f)} = 
-- -f £ (w-z) I dz a S(z-z) I dw b S(w -w)O zw Tr[U(<y)] + 



7 



7 



1 



+f e (w - z)- j dz a S(z dw b S{w - w)Tr[U(j)} 



dw b 5(w - w)O zw Tr[U(~f)} 



(CI) 



Where we have used that O wz Tr[U(j)} = O zw Tr\U(^f)}. Recall, when comparing with section VI, that we have 
omitted a factor j in the commutator of the Hamiltonians. 
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